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Abstract. Short-term variability of X-ray continuum spectra has been re- 
ported for several Active Galactic Nuclei. Significant X-ray flux variations are 
observed within time scales down to 10 3 — 10 5 seconds. We discuss short vari- 
ability time scales in the frame of the X-ray flare model, which assumes the 
release of a large hard X-ray flux above a small portion of the accretion disk. 
The resulting observed X-ray spectrum is composed of the primary radiation and 
of a reprocessed Compton reflection component that we model with numerical 
radiative transfer simulations. The incident hard X-rays of the flare will heat 
up the atmosphere of the accretion disk and hence induce thermal expansion. 
Eventually the flare source will be surrounded by an optically -thick medium, 
which should modify the observed spectra. 



1. Introduction 

In the last few years several Active Galactic Nuclei (AGN) showing short-term 
X-ray variability ha ve been reported (see for instance Ark 5 64 and Ton S180 
(lEdelson et all2002h. I Zw 1 (jGallo et all2004t ). MGC-6-30-15 (|Vaughan fc Fabianl 



2003: lPonti et al.ll2004l )h Modeling this behavior remains a complex task: owing 



to the short variability time scales observed, one has to assume small emission 
regions down to the size of a few Rschw (Schwarzschild radii). Such regions ex- 
ist within the framework of the flare model. In this model compact hard X-ray 
sources resulting from magnetic reconnection evolve above the accretion disk. 

In the following we present detailed numerical modeling of the reprocessed 
X-ray spectra induced by a magnetic flare. We investigate the influence of a 
realistic flare illumination as compared to an isotropic one and the dependency 
on the disk inclination. Using this model we attempt to derive the consequences 
for the time evolution of the observed spectra. 

Then we briefly discuss the hydrostatic balance and the optical properties 
of a the disk atmosphere below a flare. 



2. The composite spectrum of a flare-irradiated hot-spot 
2.1. Modeling 

We assume a plane-parallel model atmosphere at a typical distance of QRschw 
from the central black hole. iGollin et alJ (|2003T ) have shown that, after the on- 
set of a flare, the thermal and the ionization profile change much quicker than 
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the hydrostatic equilibrium. We therefore assume the steady density struc- 
ture of a non-ir r adiate d disk in hydrostatic equilibrium calculated according to 
iRozariska et al.l (|l999h . The incident flux follows a power law with Fi nc oc E~ a , 
a = 0.9 over 0.1 eV up to 100 keV. We choose a ratio Fi nc /Fdi s k = 144 of in- 
cident flux of the flare to the flux coming from the disk, hence the flare by far 
dominates over the disk emission. 
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Figure 1. Flare geometry: the observed spectrum depends on the incident 
angle of the irradiating photons and on the viewing angle. 0o denotes the 
half-opening angle of the flare-cone. 



In previous flare models the disk illumination was defined as semi- isotropic 

( Czernv et al.l2003 : ICollin et all2003l ) or with a constant incident angle ( Navakshin fc Kallmanl 
l2oniT >. However, the incident angle and the ionization parameter £ = 47r " F '" c of 

the flare (with Fi nc being the incident flux and nn being the hydrogen density 
at the disk surface) change significantly with the distance to the center of the 
illuminated spot (Fig^). In particular we have £ = £,o(cosQ)~ 3 if £o denotes the 
ionization parameter at the spot center. 

We have to define a maximum radius of the irradiated spot, i.e., a maximum 
half-opening angle Bo of the flare cone. We choose Oo = 60° because within the 
corresponding spot 50% of the luminosity reaching the disk is reprocessed. The 
observed reprocessed spectrum then consists of the integrated spectra of succes- 
sive concentric rings at incident angles with 0° < B < 60°. These spectra can be 
computed for different inclination angles i (measured from the disk symmetry 
axis) . 

We perform detailed radiative transfer simulations coupling the codes TI- 
TAN and NOAR. TITAN simulates radiative transfer inside o ptically thick me- 
dia (|Dumont. Abrassart. Collinl 1200(1 : iDumont et alJ 12003) determining the 
temperature and ionization structure. NOAR is a Monte-Carlo code treating 
the Compton processes. 
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Figure 2. Integrated reflection spectra for different half-opening angles Qq 
seen at a typical Seyfert-1 inclination i ~ 25°; for larger half-opening angles 
the spectra become harder. 



2.2. Results 

The reflected spectra have been computed for successive rings of the illuminated 
hot spot up to @o = 60°. While the Compton hump in the hard X-ray range is 
rather similar for all incident angles considered, the spectral slope in the soft X- 
rays rises for emission regions farther away from the spot center. The integrated 
Compton reflection for several half-opening angles seen at an inclination of about 
25°, which represents a typical Seyfert-1 inclination, are shown in Fig|2j 

Due to different distances to the flare source the reprocessing will evolve 
in time from the center of the spot toward the outer regions. This can be 
seen as successively growing half-opening angles of the flare cone. The spectra 
should become harder as the flare intensifies. The time scale of this evolution is 
determined by the height of the source above the disk. Such a developme nt was 
observed for the Narrow Line Seyfert 1 galaxy I Zw 1 ( Gallo et al.ll2003 ). 



3. Flare-induced evolution of the disk atmosphere 

After the dynamical time scale, the atmosphere settles into a new hydrostatic 
equilibrium with a total height Htot and a Thomson optical depth T es . We 
constrain the effect of a flare on the total height Htot an d the Thomson optical 
depth r es . 

We analytically recompute the new equilibrium and again include the angu- 
lar dependence of the incident flux for different distances d from the spot center. 
We find the general result that the expansion of the disk atmosphere induced 
by the flare and its optical depth depend strongly on the distance of the flare 
from the central black hole. 
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Figure 3. Top: radial dependence of the thickness (in flare height units). 
Bottom: Optical depth of the heated zone at the spot center. The different 
curves denote three values of the energy fraction f cor dissipated in the corona. 



In FigEl we show the ratio Htot/h, with h being the height of the flare 
source, and r es versus R/Rschw We choose a black hole mass of M = 10 s M , 
an accretion rate m = 0.03 (Eddington units), and a covering factor of the disk 
surface with hot-spots of $ cover = 0.1 (see details of this work in Czerny & 
Goosmann (submitted)). The different curves refer to three values of the fraction 
of flux dissipated in the corona. The tendency in all cases is that toward the 
outer regions of the disk the atmosphere reaches a higher altitude and becomes 
optically thinner. For intermediate distances the flare source eventually will be 
shrouded by an optically thick (r es > 1) disk atmosphere {Htot/h > 1) which 
should modify the observed spectrum. 
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